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Three-dimensional  Crack  Depth  Profile 
Assessment  using  Near-Field  Surface  Acoustic 

Wave  Signal  Response 


James  L.  Blackshire,  Member,  IEEE 


Abstract — A  method  for  determining  the  three-dimensional 
depth  profile  of  a  surface-breaking  crack-like  feature  is 
presented  based  on  near-field  surface  acoustic  wave  signal 
responses.  Three-dimensional  finite  element  models  were  used  to 
study  the  forward  problem,  where  the  characteristic  near-field 
scattering  of  a  surface  acoustic  wave  incident  on  a  surface¬ 
breaking  crack  was  investigated.  Experimental  validation  of  the 
modeling  predictions  was  accomplished  using  a  wedge 
transducer  for  surface  wave  generation  and  a  scanning  laser 
vibrometry  system  for  surface  wave  detection.  The  characteristic 
near-field  amplitude  response  in  reflection  and  in  transmission 
was  measured  and  modeled  for  flat-bottom,  angled,  and  curved- 
bottom  localized  notch  features,  where  a  simple  linear  inversion 
method  was  developed,  which  provided  an  effective  means  for 
characterizing  and  mapping  the  three-dimensional  depth  profile 
of  surface-breaking  crack-like  features  with  depths  in  the  micron 
to  millimeter  range. 

Index  Terms — Surface  Acoustic  Wave,  Crack  Depth  Profile, 
Near  Field  Scattering 


I.  Introduction 

N  recent  years,  the  need  for  fully  characterizing  damage  in 
aerospace  systems  has  become  increasingly  important.  The 
accurate  and  reliable  quantification  of  crack  length,  depth, 
aspect  ratio,  and  orientation,  for  example,  is  needed  for 
condition  based  maintenance  (CBM)  assessments  of  many 
legacy  and  future  aircraft  systems  [1].  Ultrasonic  sensing 
represents  a  key  NDE  method  for  2D  and  3D  crack 
characterization,  where  recent  advances  in  sensing  and 
analysis  methods  have  led  to  important  improvements  in  both 
measurement  capabilities  and  accuracies  [2-7]. 

Because  they  typically  reside  on  accessible  surfaces, 
surface-breaking  cracks  (SBCs)  represent  an  excellent  test 
case  for  understanding  and  advancing  2-D  and  3-D  crack 
characterization,  where  the  specific  features  of  a  crack  can  be 
easily  verified  [2,8].  In  addition,  the  analytic  scattering 
problem  has  also  largely  been  solved  for  surface  acoustic 
wave  interactions  with  surface-breaking  cracks,  providing  a 
sound  technical  basis  for  method  development  and  validation 
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[9-11].  Recent  advances  in  distributed  sensing,  and  point-wise 
generation/detection,  for  example,  have  led  to  improvements 
in  sensing  resolution,  where  innovative  sensing  concepts  near 
a  crack  scattering  site  can  be  exploited  [11-13]. 

In  the  present  effort,  the  use  of  forward  numerical  models, 
point-wise  laser  interferometric  detection,  and  inverse  method 
development  have  been  combined  to  provide  a  means  for 
quantifying  the  local  3-D  crack  depth  of  a  surface-breaking 
crack.  Numerical  models  were  used  to  parametrically  study 
the  effect  of  depth  variation  on  the  local  reflected  and 
transmitted  SAW  amplitude  levels  immediately  in  front  of, 
and  immediately  behind,  simulated  surface-breaking  crack 
features  with  known  depth  profiles.  A  simple  linear 
relationship  was  observed  between  the  local  amplitude  ratios 
and  notch  depth,  where  a  simple  inversion  formula  was 
developed,  and  experimental  studies  were  subsequently  used 
to  validate  the  model  predictions  and  inversion  method 
approach. 

II.  Saw  scattering  from  surface-breaking  crack 
A.  Physics  of  Scattering  Problem 

The  scattering  of  a  surface  acoustic  wave  from  a  surface¬ 
breaking  crack  has  been  studied  by  many  researchers  and  is 
well  understood  [1,14-16].  To  a  large  degree,  these  efforts 
have  focused  on  the  analysis  of  far-field  signals,  where  in 
general  the  reflection/transmission  signal  behavior  is  predicted 
to  increase/decrease,  respectively,  with  increasing  crack  depth. 
Numerous  experimental  studies  [5,17-19]  have  verified  the 
basic  features  of  these  analytic  predictions,  with  numerical 
models  [2-5]  providing  additional  insight  into  the  problem. 

The  development  of  practical  and  reliable  characterization 
methods  still  remains  a  challenge,  however,  where  the 
complex  nature  of  the  scattering  process  can  result  in 
numerous  confounding  factors  involving  the  spatial  and 
temporal  interference  of  mode-converted,  reflected,  and 
transmitted  waves  [2,18,20].  These  interference  effects  can 
cause  oscillatory  behaviors  in  the  far-field  response,  which 
can  make  detailed  crack  characterization  difficult  and 
unreliable.  In  addition,  source  variations,  transducer  coupling 
differences,  and  sensor-crack  separation  distances  have  also 
been  suggested  as  causes  for  the  response  variability  [18,20]. 
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Recent  ultrasonic  approaches  involving  measurements 
made  very  near  the  surface-breaking  crack  site  have  suggested 
that  opportunities  for  enhanced  characterization  and  reliable 
sizing  may  be  possible  [20-26].  Local  enhancements  of  the 
SAW  amplitude  near  a  surface-breaking  crack  site  have,  for 
example,  been  observed  experimentally  [13],  where  the 
enhancements  closely  follow  the  surface  crack  morphology, 
providing  a  means  for  directly  mapping  the  detailed  surface¬ 
breaking  crack  features. 

Additional  opportunities  for  obtaining  local  crack  depth 
estimates  have  also  been  suggested  [26],  where  measurement 
points  immediately  in  front  of  and  behind  a  surface-breaking, 
crack-like  feature  are  of  interest  (Figure  1).  In  particular,  the 
out-of-plane  displacement  levels  at  those  two  locations  vary  in 
a  systematic  manner  with  the  surface-breaking  feature  depth. 
For  shallow  depths  (d/ZR  «  1),  very  little  energy  will  be 
reflected  or  scattered,  with  most  of  the  energy  being 
transmitted  as  RT3.  As  depths  increase,  scattering  from  the 
front  root  (Si)  will  combine  with  the  incident  Rayleigh  wave 
(Ri),  with  reflected  energy  from  the  front  face  (Rr),  and  with 
transmitted  energy  along  the  front  face  (Rn)  to  create  an 
enhanced  displacement  field  immediately  in  front  of  the 
surface-breaking  feature.  Forward  propagating  energy 
scattering  from  the  tip  (S2)  will  reach  the  surface  immediately 
behind  the  surface-breaking  feature  to  enhance  displacements 
slightly,  before  decaying  as  depths  increase  further. 
Enhancements  of  the  energy  in  front  of  the  surface-breaking 
feature  will  continue  to  increase  as  depths  increase  until  the 
penetration  depth  of  the  SAW  begins  to  approach  the  feature 
depth  (d).  Beyond  d/XR  ~1,  reflected  and  transmitted  energy 
levels  are  expected  to  level-out,  as  incident  (Ri),  reflected 
(Rr),  transmitted  (Rti),  and  front  root  scatter  (Si)  waves 
stabilize  and  combine  to  produce  the  enhanced  displacement 
levels  immediately  in  front  of  the  feature,  while  the 
transmitted  (RT2  and  RT3)  Rayleigh  waves  and  scattered  back 
root  (S4)  energies  will  stabilize  and  combine  to  give  minimal 
displacements  immediately  behind  the  feature. 


Fig.  1.  Schematic  diagram  of  Rayleigh  wave  (R,)  incident  on  surface¬ 
breaking,  crack-like  feature  with  reflected  wave  (Rr)  from  front  face, 
scattering  from  front  root  (Si),  scattered  wave  from  tip  (S2),  scattered  wave 
from  back  root  (S4),  and  transmitted  Rayleigh  waves  (RTi,  Rt2,  Rt3)- 


B.  Numerical  Modeling  Approach  and  Studies 
A  series  of  3D-fmite  element  models  were  used  to  study  the 
behavior  of  the  Rayleigh  wave  interactions  with  surface¬ 
breaking,  crack-like  feature  as  depicted  in  Figure  2.  The 


model  was  developed  in  a  commercially  available  finite 
element  package  PZFlex,  which  is  designed  and  optimized  for 
elastic  wave  propagation  analysis.  The  model  geometry 
depicted  in  Figure  2  was  chosen  to  approximate  a  l/2-penny 
surface-breaking  crack.  The  overall  dimensions  of  the  model 
included  a  10  mm  x  10mm  cross-section  and  20mm  length, 
with  the  rounded  bottom,  ^-penny  feature  located  along  the 
top  surface  at  the  mid-line  along  the  geometry  length.  The 
feature  was  oriented  normal  to  the  top  surface  and  had  a  width 
of  25  um,  with  depths  varying  from  0  microns  on  the  left  and 
right  sides  to  a  maximum  depth  of  500  microns  at  the 
midpoint  in  the  transverse  direction. 


Fig.  2.  Schematic  of  3D  finite  element  model  geometry  (left),  and  cross- 
section  of  surface-breaking,  crack-like  feature  (right). 

The  material  properties  of  aluminum  were  used  for  the 
study  with  a  density  of  2700  kg/m3,  longitudinal  velocity  of 
6419.88  m/s,  shear  velocity  of  3039.86  m/s,  and  Rayleigh 
velocity  of  2996  m/s.  A  shear-dipole  excitation  source  was 
applied  along  the  top  left  surface  to  generate  and  propagate  a 
Rayleigh  wave  along  the  material  surface  towards  the  V2 
penny  feature.  The  separation  distance  between  the  source  and 
the  l/2-penny  feature  (10  mm)  was  chosen  to  provide  adequate 
separation  of  the  surface  skimming  longitudinal  wave,  bulk 
waves,  and  Rayleigh  wave,  which  are  all  generated  using  the 
shear  dipole  source  excitation.  Absorbing  boundary  conditions 
were  applied  to  the  all  of  the  exterior  surfaces  except  for  the 
top  surface  which  was  set  to  be  a  free  boundary.  An  example 
of  a  typical  excitation  source  is  depicted  in  Figure  3,  along 
with  its  frequency  spectrum  for  a  1 .25  MHz  excitation. 


Fig.  3.  Finite  element  1.25  MHz  excitation  signal  (top),  and  its  frequency 
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spectrum  (bottom). 

Previous  research  had  shown  characteristic  signal  amplitude 
variations  and  enhancements  which  varied  with  notch/crack 
depth  in  the  immediate  vicinity  of  the  scatter  site  [4,13,26]. 
Figure  4  depicts  model-based  signals  extracted  from  positions 
immediately  in  front  of  and  behind  a  crack-like  feature  for  a 
SAW  frequency  of  1.25  MHz  (kR  ~  2400  um)  and  for  notch 
depths  between  250  um  and  450  um  (d/^R  =  0.10  and  0.19, 
respectively).  Reflection  amplitudes  show  an  increase  of  32% 
for  a  notch  depth  increase  from  250  um  to  450  um,  while  a 
slight  decrease  in  amplitude  can  be  observed  for  the 
transmission  case. 


Omm,  NF  Reflect,  1.25MHz 


- 250um 

- 300um 

- 350um 

- 394um 

- 400um 

- 450um 


-0.0mm,  NF  Transmit,  1.25MHz 


- 250um 

- 300um 

- 350um 

- 394um 

- 400um 

- 450um 


Fig.  4.  Out-of-plane  displacements  for  1.25  MHz  SAW  frequency,  feature 
depths  between  250  um  and  450um,  and  for  position  immediately  in  front  of 
notch  (top),  and  position  immediately  behind  notch  (bottom). 

Figure  5  (top)  provides  the  characteristic  curves  for  the 
normalized  peak-to-peak  amplitude  response  versus  d/XR  for 
the  signals  depicted  in  Figure  4.  Relatively  complex  curves 
are  noticed  for  both  the  reflected  and  transmitted  cases,  with 
sharp  rising  features  below  d/^R  =  0.3,  a  peak  occurring 
between  d/XR  =  0.3-0. 5,  and  a  plateau  or  leveling  beyond  d/XR 
=  0.8  in  both  cases.  Estimating  the  depth  of  a  crack-like 
feature  using  these  curves  would  be  difficult.  However,  by 
taking  the  ratio  of  the  normalized  peak-to-peak  signals 
immediately  in  front  of  and  behind  the  crack-like  feature,  the 
bottom  curve  in  Figure  5  is  produced.  A  linear  trend  is 
observed  for  this  ratio  plot  between  the  d/?iR  =  0. 1-0.8  range, 
which  can  be  used  to  estimate  the  local  crack  depth  using  the 


equation: 


d=C,MNF„„+C2))  ,  (1) 

where  is  the  SAW  wavelength,  and  Nfrati0  =  NFT  /  NFR  is 
the  ratio  of  the  measured  peak-to-peak  amplitude  levels 
immediately  in  front  of  and  behind  the  notch,  respectively. 


d/lambda 


Fig.  5.  Normalized  peak-to-peak  amplitudes  for  out-of-plane  displacements 
obtained  in  front  of  and  behind  notch  for  increasing  d/XR  values  (top),  and 
ratio  of  the  near-field  reflection  and  near-field  transmission  plots  (bottom). 


C.  Surface-Breaking  Depth  Characterization 
The  proposed  methodology  for  estimating  the  local  depth  of 
a  crack-like  feature  uses  the  simple  linear  transform  equation 
(Equation  1),  where  Q  and  C2  are  proportionality  constants 
dependent  on  the  specific  crack  geometry,  XR  is  the  Rayleigh 
wavelength,  and  NFratio  =  NFR/NFT  is  the  ratio  of  the  measured 
displacement  amplitudes  immediately  in  front  of  (NFR)  and 
behind  (NFT)  the  feature.  Because  Rayleigh  waves  are 
inherently  2-dimensional  in  nature,  the  local  depth  in  the 
transverse  direction  can  be  mapped  if  spatially-resolved 
displacement  measurements  in  the  transverse  direction  are 
made.  Figures  8  and  9  provide  examples  of  the  predicted 
depth  profile  vs  actual  depth  profile  for  an  angled-bottom  slot 
in  aluminum,  and  a  model-based  !/2-penny  crack-like  feature 
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as  described  in  the  previous  sections,  where  good  agreement  is 
observed  for  the  actual  depths  versus  predicted  depths. 


Transverse  Position  (mm) 


Fig.  8.  Plot  of  predicted  versus  actual  notch  depth  (bottom)  for  an  angled- 
bottom  slot  using  measured  near-field  amplitude  ratios  and  Equation  (1). 


Fig.  9.  Comparison  of  model-based  Zt-penny  crack-like  feature  with 
input/actual  model  profile  depth,  and  predicted  depth  using  model-based  near¬ 
field  amplitude  ratios  and  Equation  (1). 
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III.  Conclusion 

A  method  for  determining  the  three-dimensional  crack 
depth  profile  of  surface-breaking  crack-like  features  is 
presented  based  on  near-field  surface  acoustic  wave  signal 
responses.  Three-dimensional  finite  element  models  were  used 
to  study  the  forward  problem,  where  the  characteristic  near¬ 
field  scattering  of  a  surface  acoustic  wave  incident  on  a 
surface-breaking  crack  was  investigated.  The  characteristic 
near-field  amplitude  response  in  reflection  and  in  transmission 
was  measured  and  modeled  for  flat-bottom,  angled,  and 
curved-bottom  features.  A  simple  linear  inversion  method  was 
developed,  which  provided  an  effective  means  for 
characterizing  and  mapping  the  three-dimensional  depth 
profile  of  surface-breaking  crack-like  features  with  depths  in 
the  micron  to  millimeter  range. 
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